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The Minimal Supersymmetric Standard Model has several flat directions, which can naturally be 
excited during inflation. If they have a slow (perturbative) decay, they may affect the thermalization 
of the inflaton decay products. In the present paper, we consider the system of udd and QLd flat 
directions, which breaks the U{1) x SU{2) x SU{3) symmetry completely. In the unitary gauge and 
assuming a general soft breaking mass configuration, we show that for a range of parameters, the 
background condensate of flat directions can undergo a fast non-perturbative decay, due to non- 
adiabatic evolution of the eigenstates. We flnd that both the background evolution and part of the 
decay can be described accurately by previously studied gauged toy models of flat direction decay. 
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Flat directions are generic features of supersymmetric theories. They are directions in field space along which the 
renormahzable part of the scalar potential vanishes. The Minimal Supersymmetric Standard Model (MSSM) and its 
extensions have a plethora of D and F-flat directions T|, which are lifted due to supersymmetry breaking. During 
, inflation, if their effective mass is small compared to the Hubble rate, the fields can develop large vacuum expectation 
' values (VEV) along the flat directions of the potential [1, Q . This growth is bounded above by the non-renormalizable 
(— I ' term with lowest dimension, which has the form (j/ / M'^~'^ with d > 4. If inflation is long enough, the growth will 
Q^' proceed up to (</)) - {m^M'^-^ f/^'^-'^^ If all non-renormalizable terms allowed by gauge invariance are present, 
' \ each class of flat direction will be lifted by the term with the smallest d allowed by MSSM symmetries On the 
other hand, discrete symmetries may forbid some of such terms and a non-renormalizable term with a higher d may 
(— I ' determine the VEV of the flat direction. All possible flat directions in MSSM are excited by terms with d < 9 ^. 
I— I: The formed condensate can have several cosmological implications: In the presence of phase dependent potential 
terms, the flat directions may source a finite baryon number density through the Affleck-Dine mechanism @, 0i 0- 
, It has also been suggested that they may be responsible for inflation Q . Additionally, due to the large VEV of the 
flat directions, all the fields coupled to them acquire a large effective mass, slowing down the decays they mediate 
and resulting in a small perturbative decay rate. Typically, the perturbative decay of flat direction concludes after 
00 ' ^ 10^^ rotations [l^]- These long lived flat directions also keep the gauge fields of broken symmetries (assumed to 
' be all the Standard Model ones) heavy, suppressing the scatterings among the inflaton decay products, thus delaying 
, their thermalization In addition, the energy density of the (relativistic) inflaton decay products may become 

I sub-dominant over that of (massive) flat directions. The subsequent radiation stage will then be dominated by the 
0^ : thermal distribution of flat direction decay products, rather than those of the inflaton. These effects on thermalization 
require sufhciently large initial flat direction VEVs, which can be acquired only if non-renormalizable superpotential 
^ ] terms up to d = 11 are absent (lOj . 

However, if the decay of the flat directions is controlled by non-perturbative effects^the effect on thermalization 
will be very different than the above picture. This possibility was first discussed in [12], in the framework of a toy 
^ I model based on F-term type interactions. For this model, the frequencies of the particles coupled to the flat directions 
evolve adiabatically, not allowing a resonant decay. On the other hand, it was shown in [IQ] , that the D-term potential 
provides non-trivial interactions among the perturbations through a non-diagonal and time dependent mass matrix. 
Even if the eigenvalues of this matrix evolve adiabatically, the diagonalization procedure itself may be non-adiabatic, 
due to a fast rotation of the eigenvectors. The resulting exponential decay of the condensate has a much higher rate 
than the perturbative one, giving a decay after 0(10) rotations of the flat directions. In [l^], it was also argued 
that at least two or more flat directions need to be excited for this effect to be realized. The argument is as follows: 
Since the resonant effect occurs in the D-terms, only the perturbations coupled to the VEVs through the symmetry 
generators are counted. Out of these degrees of freedom, two per broken symmetry will correspond to a Higgs and a 
Goldstone. Furthermore, two more (light) degrees of freedom will decouple, corresponding to the real and imaginary 
parts of fluctuations along each flat direction. In order to have a non-adiabatic mixing, one needs additional light 
degrees of freedom that the condensate can decay into. To formulate, the number of remaining degrees of freedom 



X 



^ The term at which all flat directions are lifted may be different for the extensions of MSSM. For example, in vMSSM, no flat direction 
survives beyond d = 6 
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present in the system will be 

/ remaining \ / d.o.f. \ 2 x ( broken \ 2 x ( 
\ degrees I I in D terms J I symmetries J I directions 

As long as this number is zero, there will not be any room for non-perturbative decay. For instance, for the typical 
cases of single flat directions, no residual degree of freedom is present On the other hand, there exist flat 

directions that are non-exclusive, i.e. they do not give a large mass to each other due to their VEVs. If the conditions 
to excite a single flat direction are present, one can expect that the whole set of flat directions non- mutually exclusive 
with that one is excited. If realized, such a case would provide the extra degrees of freedom into which the condensate 
may decay non-perturbatively. 

The longevity of single flat directions was later reiterated by the authors of , where it was also argued that for 
the non-perturbative decay of multiple flat directions, one needs some degree of tuning of the initial VEVs: Since 
different flat directions may be lifted by different non-renormalizable terms in the superpotential [l|, one may in 
general expect hierarchical VEVs. Such a case reduces to a single flat direction, which decays only perturbatively. 
The maximum amount of hierarchy that can provide a non-perturbative decay depends on the ellipticity of the orbits 
of the VEVs in their complex plane. In later works, gauged toy models with two flat directions [HI, [l^ and examples 
from MSSM [l^ were studied, each verifying that multiple flat directions may decay non-perturbatively. Additionally, 
in [1^, the fast decay was shown to be realized for a range of VEV ratios of three orders of magnitude. This range 
was found to be a consequence of the phase dependent terms introduced in the fashion of ^ . 

However, the question of whether the toy models provide a good description of MSSM flat directions needs to be 
answered. For example, in [15], the gauged toy models of four fields with only U{1) or SU{2) charges have been 
studied, yet in MSSM, no such flat direction configuration is possible and generically, for multiple flat directions, the 
field content has charges of all symmetries. Furthermore, the production of the remaining degrees of freedom ([1]) may 
be suppressed if they acquire large masses through the F-terms. Therefore, the main goal of the present work is to 
find a concrete example from MSSM for which, the decay of the flat directions proceed analogously to the gauged 
toy model case. The flat directions in the models of [l^ are decoupled at the background level, and each of the 
two VEVs evolve independently like two single flat directions. However, having independently evolving VEVs is not 
a requirement for non-perturbative decay. For instance, flat directions with coupled VEVs also have the necessary 
ingredients for decay [3, [l^ . The latter systems are much more complicated than the former ones and a precise 
answer requires extended numerical calculation. Our primary focus will be on the system of u'^u'^d'^ and QLdP flat 
directions and we will show that their decay can be described by the four field toy model of p^ . 

Additionally, we will address some issues arising from the assumption that the fluctuations along the flat directions 
are decoupled from the other modes. For instance, for the QLdF -f LLe"^ system, Ref. [Tsl] claimed that the Higgses 
and the flat direction perturbations have non-adiabatic mixings. On the other hand, for the toy models of 15], it 
was shown that the these degrees of freedom indeed decouple from the rest of the action. However, the latter result 
is a consequence of the assumption that the fields in a given flat direction have equal masses, an assumption not 
generically applicable to MSSM fields. If the fields have distinct masses, the flat direction perturbations are no longer 
decoupled from the Higgses. If these mixings are non-adiabatic, they may result in a non-perturbative decay, even 
if the counting ([T]) leaves no extra degrees of freedom. For the models we consider, we show that these mixings 
have negligible effect and the flat direction perturbations decouple as described in [lO, EB] ■ We will first generalize 
the single flat direction toy model of [iH] to have arbitrary masses and verify that the flat direction does not decay 
non-perturbatively. The approximations and methods we adopt in this simple example will provide us the necessary 
tools for the background evolution of u'^u'^cP and QLdF flat directions, which will also be studied with generic mass 
terms. 

The paper is organized as follows. In Section |TT1 we discuss different classes of multiple flat directions in MSSM over 
some examples which allow for non-perturbative decay, and determine which example is most likely to be described 
by the gauged 4- field toy model. In Section [1111 we review the formalism for the calculations of this decay. In Section 
llVi we generalize the single flat direction toy model with two complex scalar fields and U{1) gauge field, to include 
arbitrary soft masses. This provides a basis for non-degenerate mass calculation of a more complicated model, carried 
out in the next section. In Section |V] we present a complete study of u'^d'^d'^ + QLd'^ flat directions in MSSM, with 
arbitrary soft masses, where we compare the final results to the ones for the 4- field toy model in [l^. The results are 
summarized and discussed in Section IVII Finally, we include the technical steps of the calculation in the appendices 
at the end. 
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II. MULTIPLE FLAT DIRECTIONS IN MSSM 



In this section, we classify some multiple flat direction examples in the MSSM and find a case which has the 
characteristics of the gauged 4 field toy model of [T^ . A classification of multiple flat directions can be made based 
on the evolution of the VEVs. These are, in the terminology of [Tsj . i) Overlapping flat directions, where the VEVs 
are coupled to each other, resulting in a chaotic motion of the phases; ii) Independent flat directions, where the flat 
directions are decoupled at the background level and each VEV evolves independently from the others, rotating in an 
elliptical orbit in their complex plane. For example, the simultaneous excitation of LLe"^ and QLd'^ falls into the first 
category, with the VEV choice 



(M) = (r^> = |<I>|e-, {d,)^{s\)^\^e^\ (^.,) = ^|<f>|2 + e^(-+^)/2 . (2) 
On the other hand, as an example for the second class, consider the LLe"^ and u'^d'^d'^ flat directions, with VEVs, 

K> = (4) (^^e)=(M) = (0=*. (3) 

Both of these examples were studied in where it was shown that only the first case, LLe'^ and QLd'^, exhibits the 
non-adiabatic eigenvector rotation. This result can also be deduced from the counting argument ([1]): For the LLe^ and 
QLd'^ example, the field content consists of a squark doublet, a right handed squark, a right handed selectron and two 
slepton doublets, one of which is shared by the two flat directions, as seen from The perturbations of these fields 
contain a total of 32 real degrees of freedom. On the other hand, the VEV configuration breaks U{1) x SU (2) x 5f7(3) 
down to SU{2), so the remaining degrees of freedom are 32 — 2 x 9 — 2 x 2 = 10. On the other hand, for LLe'^ and 
u'^d'^d^ , there are 3 right handed squarks, a right handed selectron and two slepton doublets, so the total real degrees 
of freedom of the field perturbations is 28. The VEV configuration Q breaks U{\) x SU{2) x SU{3) completely, so 
the extra degrees of freedom are 28 — 2 x 12 — 2 x 2 = 0. Hence, for the latter case, there is no room left for the 
non-adiabatic mixing to occur. 

The previously studied multiple flat direction toy models can also be classified based on the above criteria. For 
instance, the three field model of [TB|, with D-term potential 

^ = y (l'&il'-2|<i>2p + |<i>3n' , (4) 

and VEV configuration 



$i) = |<i>|e*% {<P2) = ^ j= e^(-+-)/\ ($3) = ||.|e'-, (5) 

falls into "overlapping flat directions" class, whereas the four field model of [l5j, with D-term potential 

F = ^(|a>ip-|a>2p + |a>3p-|a>4p)' , (6) 

with VEVs, 

($l) = ($2>-$, ($3> = (*4>=$, (7) 

is a case of "independent flat directions" . 

We stress here that we expect non-perturbative decay from both systems of "overlapping" and "independent" 
VEVs. However, our goal is to find a system of "independent" flat directions to be able to use some of the numerical 
results of the toy model [lB|- The simplest such example that has extra degrees of freedom is the simultaneous 
presence of the two fiat directions, QLd'^ and u'^d'^d'^ which breaks all Standard Model symmetries^. The remaining 
40 — 2x12 — 2x2=12 degrees of freedom in the spectrum may provide the room needed for non-perturbative decay. 
Indeed, in the detailed study in Section |Vl we verified that these extra degrees participate in a non-adiabatic mixing, 
resulting in production. 



^ This is welcome, since for the delayed thermalization argument of TTl to apply, one must assume that all the gauge bosons are heavy. 
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Of course, the existence of additional light degrees of freedom is not the only requirement for a non-perturbative 
decay. As summarized in Section |T1 the ellipticity of the VEV's orbit in its complex plane determines the range of 
initial VEV ratios for which the rapid production is realized. However, the amount of ellipticity is model dependent. 
In the case of circular orbits, the non-perturbative decay of the flat directions are highly suppressed if initial VEVs 
are not comparable. On the other hand, for a model with vanishing superpotential, the motion is purely radi al 3 , [l7j 
and even a single flat direction may undergo a fast decay. For these considerations, in the toy models of [l5(, the 
choice of @] which results in an intermediate ellipticity was adopted, where quartic terms of the form X(^^^2 + c-c) 
with A cx m^/|$oP provide the angular momentum. Due to the complexity of the MSSM example, we will assume 
that an ellipticity of the same order was acquired initially and disregard such terms afterwards, since they quickly 
become sub-dominant over the mass terms once the fields start moving, and their amplitudes decrease (due to the 
expansion of the universe) . 



2uj 
2uj 



(3' = {iuj - I) a+ [ ^ - J] a, (10) 



III. FORMALISM 

In this section, we summarize the formalism of (isj and derive the conditions under which a production is expected. 
We start from the action of N coupled scalars in Minkowski space in the form 

S=^Jd^kdT] [^-'t*' - ^'t n^{T])^] , (8) 

where is an A'' dimensional vector and fi^ is a time dependent N x N real matrix. We define C to be the orthogonal 
matrix that diagonalizes fi^ through 

C'^n^C = u? (diagonal) . (9) 

It is crucial to note that in general, the matrix C is time dependent. As a consequence, the fields in the new basis 
\P = C"^^' are not the physical eigenstates of the system as long as the matrix C evolves non-adiabatically. The 
quantum evolution equations for the generalized Bogolyubov coefficients are 18] , 



with 



where, by construction, T and / are anti-symmetric, whereas J is symmetric. Additionally, the canonical commutation 
relations impose the conditions 

aa'^ - (3* 0^ = 1, a(3'^ - P* = Q. (12) 

Finally, the occupation number of z-th physical state is 

= (/3'^/3^),^^ (no sum). (13) 

The evolution equations can be seen to be constructed out of two contributions. The first part consists of the 
anti-Hermitian matrices {±iuj — /). As such a matrix has purely imaginary eigenvalues, the effect of this part is to 
rotate produced particles of physical states into each other, while preserving the total occupation number Tr(/3'*/3'^). 
The second part however is the Hermitian matrix (w'/2w — J) which causes the change in the occupation numbers, 
thus responsible for any potential production. Following p^ . we define the "adiabaticity matrix" A as 

4- ^(2^)^-4- (ri-ip). (14) 

In order to have particle production through parametric resonance, at least one component of the matrix A should 
satisfy [lB| 

|A,|>1. (15) 



1 / ^„ 1 , 1 



/, j = - V^r^i^rv^ , v = c'c', (ii) 

2 V Jlu Ju 
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The diagonal elements of this condition 



> 1 (no sum) , (16) 



requires that the eigenfrequencies evolve non-adiabatically. Notice that this is the standard condition for non-adiabatic 
evolution, valid for the case in which the produced fields are not mixed. On the other hand, the off diagonal components 



1 1 



> 1 (no sum) , (17) 



measure the non-adiabatic evolution of the eigenvectors corresponding to the physical modes. It should be noted 
that the off-diagonal components corresponding to two degenerate states correctly vanish, as the rate of change of 
the rotation between their corresponding eigenvectors is not physical. We emphasize that the condition (fT5|) is not a 
sufhcient condition for production; a detailed numerical analysis is needed to correctly determine if non-perturbative 
decay occurs. 

Throughout the paper, we exploit the large hierarchy between the TeV scale soft masses and the VEVs, by per- 
forming a series expansion in the ratio of these scales, which we denote by e. The physical modes are expected to be 
either heavy or light, with eigenfrequencies of order 0(6°) or O(e^), respectively. The scale of the rate of change for 
the background quantities are typically of order of soft masses, so differentiation of these with respect to time raises 
the order of e by one. As the matrix C which diagonalizes the frequency matrix is unitary, its leading order in the 
expansion is the e*^ term. This implies that 

r = Oie). (18) 
From pT)) . we see that Aij is of order e or higher, for 
i. rotations between two heavy modes; 

a. rotations between states that are degenerate at least at the leading order; 

Hi. changes in the frequency of a heavy mode. 

That is, for the above cases, condition (jlSp cannot be satisfied. The only components of the adiabaticity matrix that 
can be of order e° are the diagonal elements corresponding to the light modes, and off diagonal elements corresponding 
to rotations between a light mode and another mode which may either be light or heavy. In the latter case, the 
frequencies of the two modes need to be different at leading order. Terms of order in any component of T matrix 
lead to negligible A contribution which does not change the picture and will not be calculated in this work. 



IV. 2-FIELD TOY MODEL REVISITED: NON-DEGENERATE MASS CASE 

In this section, we consider the toy model with two complex scalar fields of opposite U{1) charges and study the 
effect of arbitrary soft breaking masses. Since we expect a new mixing between the light flat direction fluctuations and 
the Higgs, we wish to verify that this new coupling does not contribute to non-perturbative decay. The approximations 
and methods adopted here will also form a basis for the background calculations of a more complicated MSSM example, 
studied in Section IVl 

We start by generalizing the single flat direction potential of [3| to arbitrary masses, 

V = m?|(/.i|2 + m2|02p + Kcj^jcbj + c.c.) + ^ + q2\cj,2\^f , (19) 

where e is the coupling constant of U{1) gauge and the U{1) charges of the fields are qi = ~q2 — 1. The quartic A 

2 I 2 

term is assumed to be cx ^ , a choice compatible with the Affleck-Dine scenario. Here, $o is the value of the flat 

direction VEV when it is still frozen. The above potential reduces to the one considered in [isj in the limit mi — m2. 



^ At the quadratic level, with all the interaction disregarded, fields of equal mass are identical. 



6 



The complete action is 
with the covariant derivatives 



(20) 



= (^a^ - y q,Af,j 0, . (21) 

Here and in the remainder of the text, we use the metric in conformal time drj = R(t)dt 

ff^. = i?2diag(l, -1, -1, -1) , (22) 

where R is the scale factor. 



A. Background 



We decompose the background fields as 



(0i) = ^^e^<^-^^^) , (02> = ^^e^(^-^^\ (23) 

where we parametrized the two phases with their sum (21]) and their difference (2 AS). Since the two fields have 
opposite charge, the U{1) transformation affects only the phase difference, which we fix at AE — 0. In the degenerate 
mass case, mi = m2, the D-flatness condition — \<f>2\ allows us to have AF = 0. However, as we show below, in 
the general case, exact D-flatness cannot be attained and we need to introduce a non-zero AF. 

For convenience, we assume that the initial gauge field vanishes {{Aj) = 0) (which is kept at zero by the equations 
of motion also at later times). With these considerations, the Maxwell's equations for {A^) can be reduced to the 
constraint 

, , , 4 FAFE' 

<^°)-;f^Taf^' (24) 



whereas the background equations of motion are 

R^ 
R 



F" +[m'R^ - — - Y.'A F+-F^ cos(4 E) + Sm^R^AF 



+ (e--2Acos(4E)+ ^;ttAFT )="' 
AF" +(^m^R^-^^- E'2^ ^F"" cos(4 E) + Sm^R^F 

+ G^-2Acos(4E) + ^^¥i:^)=0, 

- ^ ^F'f -n(4E) = 0. (25) 

where we defined m? = {m\ -f m'^)/2 , 5m? = {m\ — m'\)/2 and 5m^ can also be negative. Throughout the paper, 
a prime denotes differentiation with respect to the conformal time. From the second of (|25p . we see that the D-flat 
solution AF = is allowed only if 5m^ = 0. 

Although the above equations are all we need to solve the background evolution, it is very useful to write them in 
a series approximation, to quantify the modifications due to the introduction of a non-zero mass difference 5m^ . We 
will work in the limit 

mR.T.' \ <€.F . (26) 

In other words, we expand the equations of motion ((25)l in terms of TcV/VEV ratio. For bookkeeping, we denote the 
order of expansion by e. In this fashion, we expand the background quantities as a power series in e 
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where Q can be F, E or AF. Recalling that we chose A oc m? /F"^ , its order will be e^. 

2 rj 

Determining the leading order terms of e expansion in the time derivatives is less trivial. The term l^F in the 

equation of motion for F and the term ^F AF in the equation for AF may result in oscillations with frequencies of 
order VEV. On the other hand, the rotation of the flat directions, which has a TcV scale frequency, is only sensitive 
to the average of the fast oscillations. Therefore, we carry out this averaging formally, that is, we assume that all the 
fields evolve with the small mass scale, so each time derivative increases the order of the quantity by e. This way, 
we write the equations of motion at each order of expansion. The equations (I25|) at order ©(e^) give two algebraic 
relations 



FoAF, 



4 - " — - > 4 

which are solved by AFq — 0. When this solution is plugged in, the only 0{e^) equation is 

-Fo^AFi =0, 

which removes the 0{e^) term of AF. Next, we write down 0{e^) equations. 



0, 



F^AFo^O, 



(28) 



(29) 



R" X 
-R+2 



K' + Fo 
^ (e^ Fo AF2 + 4 Sm^R^) = , 



f2cos(4Eo) 



= 0, 



^F4sin(4So) = 0, 



(30) 



where the equation for AF2 is again an algebraic one. Finally, the 0{e^) equations are 

F[' + Fi (m^R^ _ £[,2 _ _ 2 p^^'^^'^ + ^p2 [3 ^os(4 So) - 4 FqEi sin(4 So)] = , 
— F^AFa - Sm^R^Fi = , 

(F^s;)' + 2 (FoF{ - FiF^) S;, - AF3 [Fi sin(4 Sq) + 2 FoSi cos(4 Eo)] = . 
Collecting these solutions, we find 



(31) 



£2 [F^' + eF{'+0(e2)] 



[E[,' + eI]'/ + 0(e2)] 




^0 



Fo 



■F3 cos(4I]o) 



A Fo Fi cos(4 Eq) - 2 A F^ Si sin(4 Sq) - 2 SJ, S; 



|F2cos(4So) 



'R^ 



V'2 R" 

^0 ~ ir 



A 



F^ sin(4So) 




2^0 ^0 
Fo 



AFq^Si cos(4So) 



§1 (3F^Fi-FoF{)-f S'l 



^F2sin(4So)-2§SJ, 



Fo 



,(32) 



where the leading order equations match with those of the degenerate mass case [15| . On the other hand, the solution 
for AF is 



AF 



4 Sm'^ i?2 



Fi 



(33) 



This solution guarantees that at every order computed, the large mass contribution for AF cancels with the term 
Sni^R^F. To justify the assumption of neglecting the fast oscillating part of AF, we show in Figure [1] the numerical 
evolution of AF, compared with the approximate result p3p . We evolve the exact equations of motion (|25p numerically, 
giving several initial conditions for AF that are inconsistent with p3p . We find that even for initial conditions implying 
AFi 7^ 0, the evolution of AF converges to within a tenth of a rotation, after which the approximate solution is 
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valid. On the other hand, we found that the leading order terms in equations ([32)) describe the evolution correctly as 
long as AFq = 0. This latter assumption can be justified by looking at the D-term, which reads 

Vd = ^F^AF\ (34) 

Unless the leading order term in AF is at least of 0(e), the D-tcrm will be the dominant contribution in the potential. 
Conversely, for the case where the D-term is at least comparable to the mass terms in the potential, the flat direction 
is more flat than the previous case and is a more preferred configuration. For the solution ([33]) , the D-term potential 
is, 

Vo^^e^ + Oie'), (35) 
which gives a negligible contribution to the potential of the flat direction. 




FIG. 1: The evolution of AF for the first 10 rotations of the flat direction. For this plot, we evolved the exact equations 
(|25|) . and in all cases shown, the parameters are m = 10~^ e \4>in\, 5m = m/2. The curves correspond to three different initial 
conditions for AF, where ICi is A_Fin ~ 7Sm/e, IC2 is AFin = 10~^5m/e and /C3 is A^Fin = 5 (5m^/e^/|cjf>in|. We also present 
the plot of the analytic result of ([33} for comparison, which is consistent with initial conditions AF\i^ — — 2 5m^/e^/|(/!>in|. For 
the cases ICi and IC2, the leading order term in A_Fin is 0[e), which becomes sub-dominant over the second order term within 
a tenth of a rotation. 



Finally, we note that the two terms in the magnitude of the VEV, i.e. Fq and Ai^2, redshift with different powers 
of the scale factor. Specifically, once the field starts to move at 7J ~ to (or R = i?^), i^o will evolve as 



R 



3/2 



R 



(for H < to) , 



(36) 



where we denote the initial value of the VEV as 
approach. The ratio of the two constituents is 



Therefore, it is important to check the range of validity of our 



AFo 



Fn 



5m? 



e2|<^i„|2 \R 



R 



(37) 



9 



Assuming Sm ~ TeV, \(j)in\ ~ 10~^Mp and ~ 0.1, this ratio will become of order one after about 20 e-folds of 
expansion from the time when the field starts to move. Since this corresponds to 10^^ rotations of the flat direction in 
a matter dominated background (or 10^^ rotations if radiation dominates), we can safely use the expressions ([32l) and 
([55)1 in the rest of the calculation. The time at which ([57)1 becomes of order one also corresponds to the time when 
the D-tcrm becomes comparable to the mass terms, so in the antecedent evolution, the potential of the flat direction 
will be dominated by the mass terms, just like in the degenerate mass case. 



B. Perturbations 



After specifying to unitary gauge by fixing the gauge variant degrees to zero, the perturbations can be decomposed 
in terms of gauge invariant degrees of freedom as 



2R 



-S 



H + i{F-AF)a 



2R 



(38) 



where r, 6H and a are perturbations of F, AF and S, respectively. In the limit of degenerate masses, it was shown 
in [l^ that the two real degrees of freedom in Scpi + S(f>2, which correspond to the fluctuations along the flat direction 
were decoupled from the rest and did not give any contribution to non-perturbative decay. Additionally, S(j)i — S(j)2 was 
identified as the physical Higgs, which only mixed with the (also heavy) longitudinal component of the U{1) gauge 
field, and did not result in any non-adiabatic effect. On the other hand, we will see that the presence of a non-zero 
AF will cause the light modes corresponding to the flat direction perturbations to mix with the other degrees of 
freedom. Here, we will summarize the results of the detailed calculation of the action, presented in Appendix iBl 
It is convenient to decompose the vector field into transverse and longitudinal parts through 



A, 



diL , with diA, 



0. 



The two transverse vector components decouple at the linear level, with mass 

2 _ e2 (f2 + AF^) 



gauge 



(39) 



(40) 



Unlike the case with degenerate masses, the flat direction perturbations r and a do not decouple from the rest. As 
a consequence, we end up with a coupled system of r, a, 6h, L, along with the non-dynamical perturbations of Aq. 
After integrating it out and redefining the fields to have canonical kinetic terms, we end up with an action of the 
form ([8]). Disregarding the quickly suppressed terms proportional to i?" and A, we find that the eigenmasses of the 
physical modes are, in the approximation (|26p . 







KO(e'), 








ml = 


e^FS 
4i?2 


r e2 


m% = 


e^FS 


r e2 


4i?2 


^Ur^ 



Ff) Fi 



Fq Fi 



4i?2 

„2 



4i?2 



{Ff + 2FoF2 



(F2 + 2F0F2) 



3s;,2 



i?2 



(41) 



where the terms outside parenthesis, in curly parenthesis and in square brackets are, respectively, of zeroth, first 
and second order in e. We identify the heavy mode (C) as the longitudinal vector component, as its mass coincides 
with that of the transverse components ([3D]) at the order shown. The other heavy mode {D) then corresponds to 
the physical Higgs. The light modes {A, B) can be identified as the excitations along the flat direction. All four 
eigenmasses vary adiabatically and the rotation between the modes is suppressed by high powers of e. 

This concludes the study of the single flat direction toy model with two fields, generalized to have arbitrary masses. 
We have shown that the mixing between the (light) flat direction perturbations and the Higgs does not provide a 
quick rotation, thus verifying the results discussed in [To| and computed in [l^. Once the TeV/VEV expansion is 
applied, the equations for the degenerate mass case [l5| are recovered up to sub-dominant (and negligible) terms. 



V. u(f(f + QL(f FLAT DIRECTIONS 



In this section, we consider the simultaneous excitation of u^d'^d'^ and QLd'^ flat directions, which breaks 17(1) x 
SU{2) X SU{3) completely. There are much simpler examples in MSSM with less degrees of freedom and/or broken 
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Field 


Y (hypercharge) 


degree of freedom 




u" 


4 
3 


6 






2 
3 


6 




b" 


2 
3 


6 


4>i 


d" 


2 
3 


6 


(t>5 




-1 


4 


4>G 


Qc/s 


1 
3 


12 



TABLE I: Summary table of notation and hypercharges of the field content. 



symmetries, but as discussed in Section [TTl this is the simplest case of two "independent flat directions" which has 
sufficient degrees of freedom for non-perturbative decay. 

As it is apparent from the naming scheme of the flat directions, the fields that acquire VEVs in this example are a 
squark doublet, four right handed squarks and a slepton doublet. In Table |TJ we show the definitions for the fields 0^ 
as well as the choices of the specific flavors for each particle. The D-term potential for this field content is. 



1 9 1 9 

2 8^^ 



1 



3 
a—1 



1 * 



a=l 



(42) 



where gi, 52, 93 are the coupling constants, Y, aa, Xa are the generators of U{1), SU (2) and SU (3) groups respectively, 
with SU{2) and SU(^) generators represented by Pauli and Gell-Mann matrices (see e.g. [l^). On the other hand, 
the only F-term potential that gives non-zero contribution for the choices of flavors in Table |l] is 



where yd is the Yukawa coupling in the superpotential term yd Q ■ Hd d^- 
A VEV configuration which is both F and D-flat can be chosen as 

K> = (4) = (bt) = * , (rfl) - i'^e) = (.Si) - $ . 



(43) 



(44) 



where the subscripts 1,2,3 denote the SU{3) color charges. The 12 degrees of freedom for the 6 components have been 
reduced to 4 by using the 4 D-flatness conditions and 4 gauge freedoms corresponding to diagonal generators. Since 
the SU{3) charges of the components of s'^ and Q which have VEVs are different, the choice is F-flat. Furthermore, 
the total J7(l), SU{2) and SU{3) charges of the components with same VEVs cancel, so the D-term also vanishes. 
Now, we introduce soft supersymmetry breaking mass terms in the potential, 



Vsc 



(45) 



If the mass configuration has a degeneracy such that mi = 1112 — and nii — = mg, the equations of motion 
for the background is then identical to the one in [l^ for the four field toy model, with the exception of the phase 
dependent terms that provide the rotation. However, as we discussed earlier, such a mass setup implies that different 
ingredients, such as the sleptons and squarks have equal masses. In the rest of the calculation, we will assume a 
general mass configuration in the fashion of Section HVl Generalizing the configuration ((44l) . we have 



The total potential that we will use is, 



V 



1 



$4 



soft ■ 



(46) 



(47) 



where we omit the phase dependent quartic soft terms. However, we will assume that an initial angular momentum was 
already provided after the flat directions started to evolve and their source became quickly sub-dominant afterwards. 
The configuration is still F-flat. On the other hand, as we will show, the D-flatness cannot be obtained exactly 
and we can only remove 4 degrees of freedom in the background by applying gauge transformations with diagonal 
generators. Therefore, the VEV configuration has 8 degrees of freedom. In the following subsections, we study the 
background and perturbations separately and the technical details for these calculations are summarized in Appendices 
[C]and[Dl 
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A. The model and VEV configuration 

We introduce a general decomposition of the background fields in the unitary gauge as 



G + 2AG1 ,s 



$5 



V6R 
G - AGi + AG: 



G-AG1-AG2 ,s 
- e 



\/6i? ^ ' V6i? - ' - ■ (48) 

Next, we integrate out the non-dynamical temporal components of all the gauge fields. The resulting equations of 
motion are rather bulky and are discussed in detail in Appendix[C] Here, we present them in the TeV/VEV expansion. 
First, we apply the two flat direction analogue of the approximations (j26p . 

F' G' 



/F2TG2 ' ^/WTG^ 



«{i^,G} . 



(49) 



Denoting the order of this approximation by e and expanding all fields as power series in e, we find that the equations 
of motion reduce to 



£2 [F^' + eF{'+0(e2)] 

[G'^ + eG'i + 0{e^)\ 
^ [Sj; + eS'/ + 0(e2)] 

along with the algebraic relations. 



R 



R^ Sf? 



2 -^0 ^0^2 
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Fne' 



Go 
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Go m2 /J2 



2 Eq E4 



V/2 /f" 



hO(e') 



2^0 Eq^2 



Go 



-Fo -Fq ^0 
Gi , G'l , E'l 



Go Gn 



0(e2) 



(50) 



AFi = - — 

-TO 



i?2 /5m2 



Ai^2 



.9? 



l-^e + Oie)'' 



^0 5i 



AGi = — 



i?2 /,5m? 



Go V 9t 



1-^e 
Go 



AG2 = — 



Goff2' 



LrQ 



(51) 



where we defined the mass averages and differences as 

Sml = i (2 — TO2 ~ "^3 ~ 2 m| + TOg + rng) 
5m2 = 1 (2 — TO2 — TO3 — 3 + 3 mg) , 



(5m 



•'3 ' 



2 — 1 



^ (2m?-m2- 



4m| - 2m5 - 2m|) , 



,2—1 



[ml + ml + mfj , m'^ = ^ [ml + m^ + mg) , 



(52) 



The D-terms corresponding to the diagonal generators, at the leading order are then. 



y{Y) _ 


5m\ 4 


f O(e'), 


Vjj — 


5m\ 4 

2g|' 


f O(e'), 


Vjj — 


(3 (5m| - 


f (5m^)^ 4 

„2 ^ 



0(6^), 
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As in the two field toy model example, the D-flatness is approximate and of order e^, and they will dominate over the 
mass terms much later in the evolution (see the discussion at the end of Section flVAp . 



Perturbations 



We start by expanding each field as a sum of background and perturbations, 



1 



61 = ( ^ 
^3 = 7S^('^^i' a 



^^[1 + ^] 



1 + ^ 



V3 



1 + 



5d£ 
^ ^ 2 



V3 



SC2 , (5C3 



(54) 



The perturbations to the field content can be decomposed, in unitary gauge 



(5sj — —(54 + i (5i , 
<56t = -(^lo + i<57 , 



5sl = (5/i + (5/2 + (5/3 + I (551 
((5i4 - i (5i3) , 



5u% = , 
5s\ 



F2 



(Si 



i (5i 



(5ci 
(5si 



: 5fi - 2(5/5 + i<5.g2 , 
'5/4 + <5/5 + <5/6 + z 

a 



^/4 + <5/. 



((5i6 + i ^15) , 
i - (5/6 + i (5^2 . 



^2 = (56 + i ^3 , 

(5c2 = (5i7 + z(5i8 , 



Sbl = S,h+Sf2-Sh + iSgi, 
5d% = 612 + i 5q ^ 



(55) 



5S2 



i 5-2 



(5c3 

(5S3 



<5l9 + j (520 , 

-(5ii + i (58 . 



This system, including the gauge fields, has 64 degrees of freedom. From the 3-dimensional rotational symmetry of 
the Lagrangian, we know that all the transverse vector degrees (24) decouple from the remaining 40 scalar degrees. 
Furthermore, since all the gauge symmetries are broken, the transverse vector modes will acquire masses of order VEV, 
suppressing the non-perturbative effects. The quadratic action of the perturbations of the model can be decomposed 
into 9 decoupled subsystems, which are studied in detail in Appendix [D] Formally, we have. 



(56) 



Next, we diagonalize each of these subsystems and determine the physical modes. Due to the complexity of the 
model, the diagonalization procedure needs to be carried out in TeV/VEV expansion. In Table we summarize 
the results of the detailed study of the spectrum. The heavy eigenstates which have the same mass as the transverse 
vectors in the accuracy of the expansion are identified as "Longitudinal vector modes" , whereas the remaining heavy 
modes with additional terms in the sub-leading order are labeled "Higgs" . In this example, we also encounter "other 
heavy" modes, which have masses of the order of VEV, although these masses are provided by the F-terms, and 
are proportional to the Yukawa coupling yd- Additionally, we have 2 light degrees of freedom per flat direction, 
corresponding to their fluctuations. All the remaining degrees of freedom are named "Other light" . 

The first seven subsystems do not have the basic ingredients for the non-perturbative decay. However, the last two 
subsystems have rather non-trivial mixings and they indeed provide the non-adiabatic rotation of the eigenstates. 
The system S'sj,s2,d2 which consists of the perturbations of Sj, S2 and coupled to the longitudinal vector modes 



A^^''' and A'"^^' , contains a pair of light particles and three pairs of heavy particles. We find that, in the e expansion, 
the (heavy) longitudinal vectors decouple and we get 0{e°) adiabaticity matrix components corresponding to mixings 
among the light modes and the Higgses. In other words, this system has non-adiabatic rotation which may give rise 
to a non-perturbative decay*. However, to quantify this, a numerical study is required. In such a study, one still 
needs to include components of the F matrix which do not cause production, but are responsible in the rotation of 



* We note that for the system S3 j, 32,^2) the F matrix is independent of the mass differences and do not undergo any simplification even 
in the degenerate mass limit. 
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± Vector 


1 Vector 


Higgs 


Flat Dir. 


Other heavy 


Other hght 


Total 


^non diag 


16 


- 


- 


- 


- 


- 


16 


^diag 


8 


- 


- 


- 


- 


- 


8 


Sc2 


- 


- 


- 


- 


2 


- 


2 


S'c3 


- 


- 


- 


- 


- 


2 


2 




- 


2 


2 


- 


- 


- 


4 




- 


2 


2 




- 


- 


4 






4 


4 


4 






12 






2 


2 




2 


2 


8 


Sb^ , S3 > 




2 


2 






4 


8 


Total 


24 


12 


12 


4 


4 


8 


64 



TABLE II: Table summarizing the classification of degrees of freedom for the decoupled sub-systems. 



the produced quanta into other states (eg. mixings between the Higgs and the other heavy particles). This is beyond 
the scope of the present paper. 

On the other hand, the final subsystem can be simplified down to a more familiar problem in some limits, and we 
can extract a clearer picture and a more definite result out of it. This system, consisting of perturbations to 6j, S3 

and coupled to the longitudinal vector modes A^'^*-' and A^^'^\ contains two pairs of light particles and two pairs 
of heavy particles. Upon calculation of the F matrix in e expansion, the longitudinal vector modes can be decoupled 
as before, leaving a coupled system of six degrees of freedom with order e'' adiabaticity matrix elements. In other 
words, at this level, this system has more room for production than the previous system, due to the doubling of the 
light modes. The problem can be further simplified by tuning the mass parameters. In the limit of mass degeneracy 
mi = 7712 = W3 = m, 7714 = TTis = me = "in, we see that both the eigenvalues and the F matrix components produce 
two exact copies of the coupled system found in the four field toy model of [l3|, with (7I ^ 3 That is, we end up 

with two identical systems, each consisting of two light modes and a Higgs, where particle production is due to the 
mixing between the two light modes and between the Higgs and one of the light modes, as well as the non-adiabatic 
evolution of the 0{e) frequencies. A less strict tuning, introduced by having 

. ^2 {ml + ml~2ml + 2ml-ml- mj) {Fj + Glf ^ ^ ^^^^ 

gives exactly the same F matrix elements at order AAf*'. The only differences are the other combinations of mass 
differences in the sub-leading terms of the C(e°) frequencies, as well as the high order contributions to the VEVs, 
which are non-zero in the absence of complete degeneracy. These sub-leading terms however have negligible effect on 
the occupation numbers. Therefore, for the case where AAf ^ 1, even the non-degenerate system can be described 
by the results of the four field toy model. 

The function AM is time dependent and in the course of the evolution, it oscillates. It is clear that even if AAf <C 1 
is satisfied at one time, it may be violated at some other one. However, we can estimate a condition on only the 
mass parameters by approximating Fq mRFo, G'q m R G g , S n 'y mR, Sq ifiR. Furthermore, since the 
production occurs when the instantaneous VEVs are comparable [Tj, [la| , we also take Go ^ Fq in order to guarantee 
that the condition AM <C 1 is satisfied when production is expected. With these considerations, ([Ff|) reduces to a 
time independent condition involving only the mass parameters 

^Jui^ « 1 . (58) 
(m — m) 

As long as this condition applies, the results of [lH| show that the fiat directions decay within the 0(10) rotations, 
for a range of three orders of magnitude in the ratio of the initial VEVs (see Figure 4 of [IBl)- If this condition is not 
satisfied, we still expect a non-perturbative decay; however, to quantify it requires an extensive numerical evolution. 



VI. DISCUSSION 



The primary aim of this paper was to see if the extensive numerical study carried out in the framework of 4-field 
gauged toy model ^15| can accurately describe a realistic case. We have shown that this toy model, which has only 
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U{1) symmetry, two mass parameters and no F-terms, provides a very good description of an example from MSSM, 
where both u'^dPd'^ and QLdP flat directions are present simultaneously. Specifically, both cases contain "independent 
flat directions", that is, the fields in one flat direction are decoupled from the others at the background level. As a 
result, the background equations of motion for the two models have the same form. Furthermore, we found that out 
of the 64 real degrees of freedom in the problem, only the two decoupled actions Ssi,s2,d!i and Sbi,s3,d3, each consisting 
of 8 degrees, exhibit the non-adiabatic evolution of the eigenstates. At the limit of degenerate masses, that is, when 
the masses of the fields in a flat direction are identical, we found that 5'bj ^33,^3 can be decomposed into two copies 
of the coupled system in the 4 field gauged toy model. Therefore, it is safe to state that, in the degenerate mass 
case, the numerical results of the toy model provides an exact description of this part of the action. If the action 
Ssi,s2,d2 contributes to the non-perturbative decay as we expect, the resulting quanta will be different ones and their 
production will not effect the occupation numbers of generated &j, S3 and (^3 perturbations, as long as the linearized 
approximation holds. 

Another focus of the present paper was to understand the effect of different soft supersymmetry breaking masses to 
the non-perturbative production. For an exactly D-flat direction, the fields in that direction have the same soft mass. 
Conversely, if all the fields have different masses, the D-term will be non-zero and proportional to the fourth power of 
mass differences. As a reference calculation, we showed that for a single flat direction, the flat direction perturbations 
still decouple at the leading order, thus verifying that there is no non-perturbative decay in this case. The effect 
of different masses in the u'^d'^df + QLd'^ example had similar consequences for the flat direction perturbations. 
Furthermore, the four field toy model with degenerate masses still describes the subsystem S}j^^si,d^ exactly when the 
condition (|57p holds. Therefore, based on the results of [l^, we conclude that the flat directions u'^d'^d'^ + QLd"^ decay 
in 0(10) rotations, also for this case. 

One issue about the example considered here is the hierarchy between the VEVs corresponding to each flat direction. 
The first non-renormalizable operator present for u'^d'^d'^ flat direction has d = 6, whereas for QLd'^, it has = 4 [l|, 
resulting in an initial VEV ratio of ~ 10'^. This ratio, although large, may still give rise to a production if the orbits 
have enough ellipticity and the mass ratio is large enough. For instance, in the numerical results of [l^ , it was shown 
that for two flat directions with mass ratio m/m — 7.63, they decay within 20 rotations for initial VEV ratio of 10^. 
If the orbits are closer to the radial one, one might still have a non-perturbative decay with a smaller mass ratio for 
the same VEV hierarchy. On the other hand, the non-renormalizable terms are model dependent, and they may be 
forbidden by some discrete symmetries. In fact, to recover the conditions of delayed thermalization ^ll'|, one needs 
4>in ^ lO^^Mp, which requires all terms with cZ < 11 to vanish ^IQ.]. 

Although the numerical results of [lB| shows that the flat directions may decay non-perturbatively, it is still not 
clear how this effect changes the picture of thermalization. Specifically, the new quanta are produced in a resonant 
band with momenta k <m which is still non-relativistic. Since the variances are large, the gauge fields will still have 
large enough mass contributions to suppress the scatterings between inflaton decay products. On the other hand, 
tracking the evolution of the produced particle distribution is beyond the reach of the linearized calculation. One 
needs to control the back reaction effects to determine how fast the variances decrease. In such a computation, the 
distribution is likely to thermalize, possibly not in 0(10) rotations, but we expect it to be much earlier than 10^^ 
rotations that is required by the perturbative decay ^ . 

It is clear that the linear study done in the present paper, along with [l^, [l^ [l^, [l^ are limited to the stages of the 
evolution until the production is significant, and they only provide a glimpse to the beginning of the non-perturbative 
decay. For instance, we found that the u'^d'^d'^ + QLd'^ example should result in a decay within 0(10) rotations for 
the range of parameters in [i5[. However, once non-linear effects become important, the particles produced through 
different decoupled actions may interact and change the outcome of this study. Therefore, one should be cautious to 
extrapolate a non-linear study based on the toy model to a realistic one. 

The logical step to be taken next is to include higher order terms and study the effect on the decay time and 
thermalization of the produced quanta. In a recent study the non-gauged toy model of [l^l was evolved on a 
lattice using the ClusterEasy code [2l[, verifying that the non-perturbative decay is still realized within the first few 
rotations. Another interesting result of was the calculation of the gravity waves, sourced by the quick decay of 
the flat directions. Their spectrum was found to fall naturally into Hz-kHz range and depending on the initial VEV 
of the flat directions, may potentially be within the reach of upcoming experiments, such as Advanced LIGO. For 
more realistic models with gauge flelds, the resulting spectrum may be different, but since the mass scale of the flat 
directions is of order TeV, it will have a frequency range similar to that in the toy model. In a future work, we will 



^ Indeed, our preliminary results from lattice simulations of non-gauged toy model |1C)I| up to 0(100) rotations show the beginning of 
thermalization of the produced particles, i.e. the initial distribution with momenta fc < m starts to extend toward high momentum 
region. 
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address the effect of back- reaction and gravity wave production in the framework of gauged models. 

However, there are still some problems left that can be dealt with analytical tools. In it was shown that 
the non-perturbative decay of flat directions may have an observable effect through the amplification of curvature 
perturbations. In the context of natural supergravity inflation [235, Ft-^f- showed that the non-perturbative decay 
of flat directions allows the inflaton preheating to be realized in these models. On the other hand, even at the 
linearized level, we do not have a complete numerical study of the decay for "overlapping flat directions" , although 
we have approximate calculations showing the non-adiabatic rotation occurs [iGj . Specifically, the toy model of 
[l^ is qualitatively different than the ones considered in the present work, as well as fl5|. in the sense that the flat 
directions are coupled at the background level. A numerical analysis in the fashion of [T^l would be very useful in 
understanding the time scale of the decay and the range of hierarchy that allows production. There are many examples 
with "overlapping flat directions" in MSSM which provide the necessary ingredients for non-adiabatic evolution (e.g. 
[l^), and it is an interesting challenge to find simple toy models that correctly describe at least parts of these examples. 
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APPENDIX A: SIMPLIFYING A COUPLED SYSTEM 



In this section, we summarize the steps for transforming any given action to the form ([5]) and conclude by discussing 
the order in the TeV/VEV expansion (|49)) that we need to calculate. We start from an action of iV coupled scalars 

[ d^kdrjCk: (Al) 



where the Lagrangian density in Fourier space Ck is, generically 

£k^^ Ti 0' + -I- 0t ^, _ ^2 0) . (A2) 

Ti and flf are symmetric and time dependent N x N matrices, whereas the real matrix Ki is in general asymmetric. 
Due to the isotropy of the problem, these matrices are invariant under the reflection operation k —k. We now 
proceed to simplify the form of this action to that of coupled oscillators in Minkowski background. We first diagonalize 
the kinetic terms with orthogonal matrix Ri 

Rj Ti Ri = Ta (diagonal) . (A3) 
We then write the Lagrangian in the new basis ip = Rjtj) a-s 

£fc = i T2 ^' + V''^ K2^ + ^^ ti^' - V'^ ^) , (A4) 

with 

K2 = Ri Ki Ri + Rl Ti R[ , ^l = Ri Ri - R'l Ti R[ - 2 R'^ Ki Ri . (A5) 
Next, we rescale the fields such that the kinetic matrix becomes unity 

R2 = {T2)-^'^ = R^ . (A6) 
In the rescaled basis x = ^2^^^: the Lagrangian reads 

Ck = I ix'h' + X'^ Ksx + X^ Kj x' - X^ X) , (A7) 

with 

if 3 = i?2 i?2 + i?2 i?^ , nl = R2nlR2-R'2T2R'2-2R'^K2R2. (A8) 
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By adding to the Lagrangian, the boundary term 

^^^^^-\^[xHK, + Ki)x], (A9) 
we obtain a more convenient form, where the matrix which mixes the fields to their derivatives is anti-symmetric 

r-k^^ix'h' + X'^ K^X-X^ K^X' -X^^lx) , (AlO) 

where 

K, = -Kl ^\{K,^ Ki) , nl^nl + l {K, + Ki)' . (All) 

Finally, to remove the mixing matrix K4, we do a further transformation with 5* = i?5 satisfying RJR'^ = K4. 
With this rotation, we obtain an action resembling coupled oscillators in Minkowski space ([8]) 

S=^Jd^kdri [*'^«'' - {nl + KIK4) Rm . (A12) 

As outlined in [15], the construction of the matrix T — C"^ C that is needed for the computation of evolution equations, 
does not require the explicit knowledge of i?5. Diagonalizing simply the matrix (r2| + KJK4^^ through 

{VlI + KIK4) ^ = Lu^ (diagonal) , (A13) 

one finds the relation C — R^S^ which gives 

r = e^e' + e^^4?. (A14) 

For the example in Section |Vl the action is very complicated and even the first steps where the kinetic matrix Ti 
is diagonalized cannot be done exactly. Therefore, for some of the subsystems we encounter, we need to apply the 
expansion (|49|) at the level of eq (jA2|. In Section ITTTl we have shown that the matrix F is needed at the order e. 
Going backwards from (jA14p . we see that to calculate this matrix at the required order, the matrix ^ is needed at 
order e°, all the Ki matrices at order e, and matrices i?2, ^i, fi? and are needed at order e^. 



APPENDIX B: CALCULATIONS FOR 2 FIELD MODEL 



In this section, we summarize the calculation for the perturbations in 2 field model of Section [IV Bl The quadratic 
action for the transverse vector modes is immediately decoupled: 



The remaining action consists of a system of four coupled real fields, 



+ {d^Aof - 2 {d,L') (5,Ao) + -^{F^ + AF') Al 



~2nlsHrSH +4\{F^ - AF^) sin(4 ^) {F r ~ AF Sh) cr 
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F^ + AF2 
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{AF r ~ F5H)Ao + 2eF AF {d,L) (9,cr) 
[Fr- AFSH)-2eFAFAo a') 



(Bl) 



(B2) 
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where, 
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Contrary to the degenerate mass case {6m ~ 0, AF = 0), r and a, which were the perturbations of flat directions, 
are now coupled to the Higgs and the longitudinal vector through the non-zero VEV difference — {412) ■ 

We then expand the fields in terms of plain waves, and solve the constraint equation for the non-dynamical degree 
Aq, which yields 



^0=1 + 



(f2 + AF2 
4P 



, eFAF , e (F2 - AF2) E' 
L' + — - — a' - .\ (AFr-FSH) 



A:2 



/c2 (i^2 _,_ AF2) 



(B4) 



Using the expression above in (IB2p . the action becomes of the form (jA2p . By a series of redefinitions as described 
in Appendix \^ we obtain the form ([5]) . The explicit expressions for the matrices K4 and flf are too involved for 
presentation. However, one can write them as an expansion series in using the e expansion (|26p and background 
expressions ([5^ and ([55)1 . The non-zero components of K4 at 0{e) are 



[K 



4J13 



[K 



4J14 



V2 
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with non-zero fl^ components up to 0{e^ 



[^4] 11 



4J13 



34 



e + m^R^-j:'^^-^ + ^ FS cos(4 So) 



[a 



4Jl4 



3AF2 sin(4So) 



2V2 



2 n2 



[^4]22 



A;^ m-^ i?^ 3 S^,-^ -f — (F{= + 2 Fq F2 j 



= ~-Fo~{^FoF, 
1 



:|^-^F^cos(4I]o) 



3A 



- Fi^ + 2 Fo F2 - — - — Fo^ cos(4 



8 



2i? 



R" 



3A 



_ i?2 „ _ j.,2 _ _ + 2 Fo F2) - — - — F^ cos(4 So) 



(B6) 



In the above expressions and elsewhere, terms outside parenthesis, within curly brackets and square brackets are, 
respectively, of zeroth, first and second order in e. 

To determine the adiabaticity conditions, the only remaining step is to diagonalize the fl^ + Kf K4 matrix. The 
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eigenfrequencies of the physical modes are 



+ R 



2 d2 



2 d2 



+ R 



R 

Rr_ 

R 



■'D 



■ Fq Fi 



■ Fq Fi 



k^ + -{F^ + 2FoF2) 



^2 + _ (Ff + 2 Fo F2) + i?2 + 3 



(B7) 



In the first two lines, we have kept the sub-leading terms proportional to A and R" for comparison with the degenerate 
mass case. At the given order of expansion, the two modes coincide with the perturbations to the flat direction in 
[l5| |. The last two lines show the frequencies corresponding to the longitudinal component of the U{1) gauge field and 



the physical Higgs, respectively. The transformation matrix ^ which diagonalizes flf 



matrix is 



V2 



( V2sin(2Eo) -V2cos(2I]o) \ 

Q ^/2 

cos(2I]o) sin(2So) -1 

Vcos(2So) sin(2So) 1 



0{e) 
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Using this, we compute the T matrix at order e, 



r = e 



/ s;, \ 

-Sf) 



V / 
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Therefore, the only mixing is between the two light modes. However, shortly after the VEV starts oscillating, the 
two light modes become degenerate due to the suppression of the term proportional to A. Based on the arguments in 
Section Hm the leading terms in the adiabaticity matrix A (fT4)) are then at least of order e. Just like in the degenerate 
mass case, this model has no non-perturbative decay channel. Additionally, due to lack of mixing in the T matrix 
between the flat direction perturbations and other modes, we conclude that the flat direction perturbations decouple 
from the rest of the system at the given order in e expansion. 



APPENDIX C: BACKGROUND CALCULATIONS FOR THE u'=d'=d'= + QLd" EXAMPLE 

In this appendix, we summarize the calculations for the background quantities in the realistic 2 flat direction 
example in Section IV Al For clearer notation, we adopt the following decomposition: 

<I>i = ^e^^S $2 = ^e^^^ ci>3 = ^e'^% 
V6i? V6i? Vqr 

$,^_^e^^. <j, G^g.s. ci,,^_^e'^% (CI) 
V6i? V6i? V6i? 



19 



Without fixing the gauge freedom, we write the background action as 



R 



1296 

A 

432 

9l 



144 



Rf^ 
R 



T.? G 



{Gl - Gl 

2 



12 



+ 



(Ff + Fi-2 Fi + Gl-Giy + Z {Fi - Fj + G? - Gl)' 

4^^' (16 Fl + 4 F| + 4 + 4 G? + 9 G^ + G^) + ll (A^^^^ + Ap^' + A^^^') {Gl + Gl) 
(4^1)^ + 4^^)' + (f 2 + Fi + Gl+ Gl) + (4"^)' + 4"=)') {F^ + Fi + Gl + Gl) 

(A^ + 4^^)') (F| + Fi) + (Ff + Fi + 4F| + G? + Gl) 



+ j=AyA^^^' (Fi-Fi + Gi + Gl) 
9192 .(^3) AY) 



18 U^)4^) (3 Gl + Gl) - f^A^^^ (V3 4'^^ + 4'-') Gl 
+ ^^^(^) [734^^^ (4 Fi + 2 F| - 2 G? + Gl) + 4^^^ (4 Fi - 2 F| + 4 - 2 G? + Gl) 



+ 



^ (51 + V3 53 4'^0 ^3 ^3 - ^ (2 91 4^^ - 3 9s 4'^^ - V3 53 4'^0 ^1 

y -2 -2 - ^ (51 ^''^ - 352 + 353 4'^^ + V353 4'^0 ^3 . (C2) 



i(5i4'''-52A^))G^S^ 



where ^41^^^ is the VEV of the gauge field corresponding to generator G. Integrating out the non-dynamical Aq fields, 
the vector fields corresponding to non-diagonal operators are constrained to vanish. The ones of diagonal components 
are 



l(A3) 



1 

91 
1 

92 
1 

93 



2FiFi 



FiiF, 



FiFi + Fi{Fi + Fi) 

2FiFi^Fi{Fi+Fi) ^, 
FiFi + FiiFi+Fi) 



G^Gg + Gj {Gl + Gl) 
3G? (Gi-Gi) 



ri2r<2 
^2^3 



Gl{Gl + Glf 



FiFi + Fi{Fi-2Fi) ^, 2GlGl-^G\ (Gj + Gj) - , _ 
FiFl + [Fi + F|) + G^Gi + G\ {Gl + Gl) ^ """^^ 



l(A8) _ 



\/353 



F|F|-Fj^(5F|-4f|) ^, + 2GiGi-Gj(Gi + Gi) ^. _ 



F|F| + Ff(F| + F|) 



GiGi + Gf(Gi + Gi) 



(C3) 



where, 



S = ^ (El + S2 + S3) , S = ^ (Si + S2 + S3) , 
ASi = i (2 El - S2 - S3 - 2 El E2 -h S3) , AS2 = ^ (2 El - S2 - E3 - 3 E2 3 E3) , 
AE3 = ^ (Ei - 2 S2 + S3 + 2 El - S2 - E3) , AS4 = ^ (Si + 4 E2 - 5 S3 -h 2 El - S2 - S3) 



(C4) 
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We replace the vector fields in (|C2[) using the above constraints, and obtain the following background equations: 

/ r>ll \ 2 

F{' + imlR^ - — \Fi + ^Fi{'iFl-2Fl-2Fl~2Gl + ZGl^Gl) 

(? q F^ 
+ § F, (2Ff - f2 - Fl ^2G\-2GV) ^ = , 



F^' + Ln\R^-?^\F^~^^F^{^Fl~2Fl-2Fl~2G\^'iGl-G, 



„2 Q p4 p p4 

36 ^ ^ 2 3 1 3; j^2^2 + ^2 (^2 + ^2^l2 



3 - 



F^' + (m^i?^-:^jF3-^F3(4Fi^-2F^-2F^-2G^ + 3G^-G^ 



^ F3 (-Ff - F| + 2 Fl - G? + G2) "^nnP^ ^2 ^ 



G'/ + [m\B'-?!^G^-^^G^ (4F2 - 2 F| - 2 F| - 2 G? + 3 G^ - G3; 



^ Gi (2 f2 - F| - Fl + 2 G? - 2 G^) ^^^^^^ .S'^ = , 

36 ^ ' ' [GiGi + G? (Gi + Gi)]' 



G^' + (mii?2_:|ljG2 + ||G2(4F2-2F|-2F|-2G? + 3G2-G2) 
G2 (Gi - Gi) 9GfG2G| ^ ~ ,2 ^ 



[GiGi + G? (Gi + Gi)] 



G'3' + [jnlR^-^^G^-$-G^ (4F2 - 2F| - 2F| - 2G? + 3G2 - G^) 



if G3 (Gi - Gi) + ll G3 (-2Ff + Fl + F32 - 2G? + 2Gi) 9GtG|G3 ^ 

24 2 3; 3g 3V 1 2 3 1 37 ^^2^2 + ^2 (^2 + G2)]2 

-pi T?1 T?1 \ ' / /^2 /^2 /^2 



.J^Ii^l + ^^iM^^I + ^D^'J UiGi + Gf(Gl + Gi)^'j °' 

where the equations of motion are independent of the gauge variant combinations AS^. Finally, with the definitions 

Fl = F + 2 AFi , F2 = F - AFi + AF2 , F3 = F - AFi - AF2 , 

Gi = G + 2 AGi , G2 = G - AGi + AG2 , G3 = G - AGi - AG2 , (C6) 

and applying the expansion (|49p. along with the assumptions AF^'/AF^ — ©(e), /S.G\/ /S.Gi = 0(e) (see the discussion 
in Section [IV Ap . the equations of motion (|C5[) give ([50)1 and ([?T|) . 



APPENDIX D: QUADRATIC ACTION FOR THE + QLd" EXAMPLE 

In this appendix, we summarize the calculations of the quadratic action for the u'^dPd'^ + QLdP perturbations of 
Section IV Bl After fixing the unitary gauge, there are in total, 64 degrees of freedom, as well as 12 non-dynamical 
degrees. As mentioned in the main text, this system is very complicated, but it is possible to pick out decoupled 
subsystems. We study the spectra and adiabaticity conditions for each of these separately in the subsections below. 
In this following, we immediately integrate out the VEV of the non-dynamical degrees (^q'^-'), so that A^j^"^ denotes 
the perturbations to the gauge field corresponding to generator G. 



1. Subsystem 5'^°'^ diag)_ -j^ansverse vectors - non-diagonal generators 

One of the immediate sub-systems of the complete quadratic action is the part containing the transverse vector 
degrees corresponding to non diagonal generators. These decouple from the rest, as well as from each other, each with 
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a very similar action differing only in the mass terms: 



1 



A 



(G) 



(G)^(G) 



(Dl) 



where all the masses are of order VEV, 
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2 _ 93 



12 

2 

(Ay) - 



"^(A4) = "^(A5) 



12 



_ .92 /'^2 

"12 



Gl) 



Gl + G\ 



3) ' 



(D2) 



This part of the action decouples 16 dynamical degrees of freedom from the rest. 



2. Subsystem S^'^^': Transverse vectors - diagonal generators 



We now move on to the transverse part of the vectors corresponding to the diagonal generators. Although these 
are decoupled from the rest, they are coupled to each other. The action can be written in the form 



(D3) 



where 



/((f^s) /l(A3) l(A8 



(D4) 



The mixing of the modes are due to the non-diagonal mass matrix components 



(M^)^^ = 4(l6Ff + 4F| 
(A^^)i2 = "^(3Gi + G^ 



4F2 



1) 

(^')i3 = "5f (4^1 +2^1 

— 3193 

2 



2G2 

-h4F. 



-4G2 



-9Gi 



/'/l,f2^ _ 9233/^2 



-2G? + Gi) 



VJ24 
iM^)3, 



12 V3 
■2 

12 I" 1 

^ (Ff - Fi + Gl 



G\ 



(^V)22 



- T2 l"-^: 



Vj44 



12 \/3 
36 V^l 



-Gl) 



4F| 



G? 
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However, as discussed in Section IIIIl this mixing does not result in a non-perturbative production of quanta. All 
the gauge symmetries are broken, so the eigenvalues of the mass matrix are all of order 1 in an e expansion. A 
straightforward way to verify this is the calculation of the determinant of My at order e°, which turns out to be 
non-zero. 

This part of the action decouples 8 more dynamical degrees of freedom. 

3. Subsystems Sea and : Scalar modes decoupled from the gauge fields 

The part of the action containing perturbations of C2 is 



(D6) 



where 



^2 F?Fi + [Ff + Fi)Fl^, ^ GlGl-2[Gl+Gl)Gl ^, 



FiFi+FHFi + Fi) GlGl + GUGl+Gl) 

AFI + 2 F^ + 2F^ + 2GI-3G 



2 

2 r/2 , Vd TTi2 



92 . 



216 

9l 
36 



g {Gl - + ^ {Fl -2Fi + F| + G? - Gl) - 



G 



(D7) 
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With these, the frequency matrix n| defined in Section IIIII is diagonal and degenerate with eigenvalues = fc^ + 
+ K^^. For non-zero yd, these modes are heavy, so the components of adiabaticity matrix A (eq. [T^]) are at least 
of order e. Therefore, non-perturbative effects are suppressed for this system. 

The second such system that does not couple to the gauge fields consists of the perturbations to C3 , with an identical 
action with {Sn, Sis) ~^ i^ig, S20) , Vd and F2 F^. The main difference of this system from the previous one is 
the lack of C(e°) terms in the frequency, that is, the eigenmodes of this system are light, with 



Sni2 



3Sr 



(D8) 



Again, due to the degeneracy of the states, the condition (flSl) cannot be satisfied. Also, the eigenfrcqucncies evolve 
adiabatically, so we conclude that there is no non-perturbative production from these systems. 
These two parts decouple four more dynamical degrees of freedom. 



4. Subsystem Se.,cx- Perturbations to e and c\ 



This subsystem consists of perturbations to e and ci, coupled to the longitudinal vector degrees of non-diagonal 
generators of SU{2). The action, in Fourier space, is 



(i4-)p + i 



A 



('T2)|2 



i|G§) (G3G2 — G2G3 
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{Gl 



GD 
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gl {3Gl - Gl) (4F2 ^2Fi - 2F| - 2G? + 3Gi - G; 



216 (Gl + Gl) 

gl G\ + 6 GlGl + Gl ^ gl Gl (-2 Fl + Fl + F| - 2 Gf + 2 G 
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Gl 
'-2FIFI 



Gl 
Fl[Fl 



nn 



F?{Fi 
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Gl + Gl 



3G? (g: 



Gl) 



GlGl + G\{Gl + Gl 



3 G\G2G^ 



\G2G5 + Gj (G2 + G3) 

+ (M5-«6+C.C.) 
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-2F|F| 



Fi (K 



FM + FliFi + Fi) 



3 (^2 — ^3) 
GlGl + Gl {Gl + Gl 



+ 
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2 73 ^Gl + G 



\^^o "16 



4(^2) c* 
■ri-n Oik 



[-2FiFi + Fi{Fi+Fi)] {Gl-Gl 



\ FiFi + Fi{Fi + Fi) 

92 {Gl - Gl) (G2G'2 + G3G',) 



3G? (g: 



Gl) 



G'iGl + Gi (Gl + G'i) 



2V3(Gi + Gi)^/' 



Next, we integrate out the non-dynamical degrees ^q'^^'' and A^q'^\ then apply the prescription of Appendix 1X1 by 
expanding in e series, to finally get the form (|A12p . In the e expansion, the mixing matrix K4 is 0{e^), and the 
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frequency matrix il^ is in block diagonal form, 



(DIO) 



where Sl^ is a 2 x 2 matrix. In other words, this action with four degrees can be separated to two identical systems 
at this approximation order. The components of the frequency matrix are 
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The eigenfrequencies of this system consists of a pair of 

52 / 92 GqGi 
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92 '-^O 
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The pairs of frequencies (A) correspond to the longitudinal vectors, as they coincide with the mass of Ap' and 
(jPip at the given order. The remaining pair of modes (B) are the Higgses. As all four of the degrees are heavy, we 
conclude that this system does not contribute to non-perturbative decay. 
This part decouples four more dynamical degrees from the rest. 



5. Subsystem 5*6, 



Perturbations to and sg 



This system is very similar to the previous one and consists of the perturbations to 62 ^-nd S3, along with the 
longitudinal components of vector fields corresponding to S'C/(3) generators Ag and Ay. The action in Fourier space 
is, 



Sh,s, = IJ d'kdJ\S 
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After integrating out the non-dynamical degrees and Aq ''' , we repeat the steps in Appendix \X\ and obtain the 

action of form ([5]). As in the previous subsystem, the matrix K4 is of order 0{e^) and ill is of form (jPlOp . with 
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The eigenfrequencies consist of two copies of 
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Again, lj\ coincides with the e expansion of the frequency of aI^""* and A^^^''^ (jPip . so it corresponds to the two 
longitudinal components, whereas lo^ are the frequencies of the two Higgses. Neither of these heavy modes will 
contribute to non-perturbative production. 

This system decouples another four dynamical degrees from the rest. 



6. Subsystem Ssf,Sg'- Perturbations to VEVs and longitudinal vectors - diagonal generators 



This system consists of the perturbations to the field components with non-zero VEVs, coupled to the longitudinal 



A 



(A3) 



vectors corresponding to the diagonal generators, namely, dfi (6 degrees), Sgi (2 degrees) and A^^\ 
^(As) rpj^g initial action is very long for presentation. Furthermore, even the zero order terms in e approximation are 
too bulky, so we will describe the spectrum and comment on how we identify the flat direction excitations. 

After integrating out the four non dynamical degrees, we need to expand the matrices in (jA2j as a series in e, and 
repeat the steps in Appendix \K\ to get the form (jA12[) . In the end, we find that at the relevant expansion order, two 
systems of each 4 degrees can be decoupled from the rest. Similar to the check we performed in Appendix ID 21 we 
calculate the determinants of their QI + KJK^ matrix at order e°, which shows that all eight of these modes are heavy, 
i.e. the system contains four Higgses and four longitudinal vector modes corresponding to the diagonal generators. 
What remain are two pairs of light modes, with leading order eigenfrequencies 
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This system is actually the analogue of the coupled system for the 2 field toy model in Section |B] As in that case, the 
perturbations to the field components with VEVs are coupled to the longitudinal gauge fields of diagonal generators. 
The only light modes in this action are then clearly the perturbations of the two flat directions. Indeed, in the 
degenerate limit Srrii = 0, -Fq = F , Gq = G, Sq = S, Sq = S, the action containing the combinations 



F F 
r = — (cos E (5/i — sin S Sgi) , a = — (sin E Sfi + cos S Sgi) , 



r = — I 
2 V 



cos S 5/4 — sin S 5g2 



2 
G 

(7 = — ( sin S (5/4 + cos S 5g2 
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immediately decouples from the rest and one recovers the action described above at the leading order in e expansion. 

To summarize, this subsystem consists of 8 heavy fields which decouple from the rest, and do not give rise to non- 
perturbative decay as described in Section lTTTl The remaining light modes correspond to flat direction excitations, with 
adiabatically evolving frequencies. Although there is a non-zero T matrix, the light modes that mix are degenerate, 
so these do not contribute to production either. 

This part of the action eliminates 12 more of the dynamical degrees of freedom. 



7. Subsystem 5'sj,s2,d2" Perturbations to sj, S2 and 



This subsystem consists of the perturbations to the field components sj, S2 and 1^2 j that is, 5i, i € [1,6], along with 
the longitudinal components of the vector fields corresponding to S'C/(3) generators Ai and A2. The decoupled action 
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93 G3 



(4'^'<5^ + 4''^^S', - {d,A['^'>m6,) {d,A['^^mS,)) - mi {Si + SI) 



%/3 

+^ {4'^^S', + 4'^^S', - {d,Ai'^^M5,) - {d^Ai'^^MSe)) - ml {Sj + Sl) 
^ ial-^vD p^Q^ ((5i(52 + SiS^) + ^F^G, {S.Ss + S^Se) - ^dGg {S2S3 + S^Se) 
+2 X.I (^1^4 - W + 2 Xs, {S2S'^ - S'^Ss) + 2 Xd, {SaS'e - S'^Se) | , (D18) 

with the mass parameters, 

ml ^ ^{F^ + Fi + Gj+Gl), 

ml ^ mlR'' + y^Gl + ^{-^Fi + 2Fi + 2F^ + 2Gi-SGl+Gl) 

+ 1 {2Fl + 2Fi-Fi + 2Gl-2Gl)-xl-^, 
ml ^ mlB? + I F| + ^ {-AF^ + 2F| + 2F| + 2Gf - + Gl) 

+ f| {-Gl + Gl) + f^{F!-2Fi + Fi + Gl + 2GD-xl-^, 
2 = m^i?^ + ^ {-AFl + 2F^ + 2F^ + 2G\-ZGl + Gl) 

+ 1 {-Fl + 2 F| - Fl + 2 + G^) - x^, - § , (D19) 



m 



and the couphngs, 
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We then expand the fields in terms of plane waves and integrate out the non-dynamical degrees Aq'^^-' and Aq^^\ 
recovering the form (|A2[) . As in the previous case, we apply the prescription in Appendix \X\ using e expansion from 
the start. The resulting + KJK4 matrix is of the form 



A B 

-B A 
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where A and B are 4x4 real matrices with symmetries A^ = A and B^ = —B. The eigenvalues of such a matrix 
consist of two copies of the eigenvalues of the Hermitian matrix A — iB. A detailed discussion on diagonalization of 
this matrix form is given in Appendix |E1 The eigenstates of this coupled system contains a pair of light modes, 
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a pair of heavy modes with order term proportional to the Yukawa coupling. 
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a pair of Higgses, 
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and finally, a pair of longitudinal vector components, 



9l 



f (^0^ + Gl) 



k' + § {Ff 
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coinciding with the masses of transverse vectors A^^^-* and from (jPip at the given order in the expansion. The 

two eigenvectors corresponding to each distinct eigenvalue at order e*^ are of the form. 



{j=A,B,C,D) 
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where aj parameters are arbitrary rotations between the eigenvectors of a degenerate eigenvalue pair. Although there 
is an additional degeneracy at zero order between heavy modes (G) and (Z?), the bottom blocks of the eigenvalues 
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turn out to be proportional to the upper blocks. The vectors are, at order e°, 

Fo 



■WA 



wc 



(0,1 ,fA(.Fo,Go) ,/a(«J^o,«Go)) , 
^^-^(0 ,1 ,/b(Fo,Go) , fB{iFo,iGo)) . 
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Go 



WD =fl, 0,0,0), 



(0,1, fc{Fo,Go) , fc{^Fo,iGo) ) 
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where the functions are defined as 

fAiFo,Go) 

/b(^o,Go) = 
fc{Fo,Go) = 



^ Fi + ^F,^ + 8Gt / -Fo^ + 2Gg + v/F4 + 8Gg 
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V2 Fo + 2 G2 + VF4 + 8G4) /F^Ts^ 
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Using the matrix at order e and the above eigenvectors, we calculate 0{e) terms of F matrix. The non-zero 
components are 



ri2 

Tie 



-a^ -I 
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-r25 = 
— 

-r45 = 
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We immediately see that the longitudinal vector modes (7, 8) are decoupled from the rest of the system and have only 
a mixing term between themselves, which does not contribute to production due their 0{VEV) mass and degeneracy. 
Since the physical quantities are unaffected by the choice of (see Appendix |E] ) , it is useful to study the T matrix 
for a suitable choice. As our main focus is to check for non-perturbative production, we look for a choice of these 
parameters which removes the non-adiabatic mixing of the eigenstates completely, that is, that makes the adiabaticity 
matrix A (eq fM)) to be of order e at least. From the discussion in Section Hill the only non-diagonal components of the 
adiabaticity matrix which has non-zero 0(6") terms will be the ones which involve the mixing of the light mode (1, 2) 
to the Higgs (5, 6). If, in this system, the non-adiabatic rotation of the eigenstates is a spurious effect, we should be 
able to remove the Fis and Fig (or, equivalently, F25 and r26) components by proper choice of rotation parameters. 
However, from (|D29[) . we see that it is not possible to make both of these components zero simultaneously. Therefore, 
we conclude that the non-adiabatic mixing of the light modes (1,2) to the heavy modes (5,6) is a physical effect, 
which cannot be removed by exploiting the freedom in the eigenvectors. In addition, the light modes' eigenfrequency 
evolves non-adiabatically, so the diagonal condition ([16]) may also contribute to production. 
This subsystem decouples a further 8 degrees of freedom. 
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8. Subsystem St^, 



Perturbations to sj, S2 and 



This subsystem consists of the perturbations Si, i ^ [7, 12], along with the longitudinal components of the vector 
fields corresponding to SU{3) generators A4 and A5. It is very similar to the one discussed in the previous subsection. 
In fact, the action for system Sbj^^g^^d^ can be obtained by 



&I, S3, d; 



= s 



n,s,,d, [S, ^ 4-^') ^ 4-^'+^), F2 ^ F3, 7712 ^ m3, yd ^ o) 
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The calculations for the action proceed the same way as the previous case, up to the point where we have the action 
of the form (|A12p . When e expansion is applied, H,^ matrices can be related by 



SI, S2, 



[dml -dml, ya O) 
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whereas the K matrix, which is independent of Sm^ and yd stays the same. 

Although we have the simple relation between this system and the previous one, the solution of the eigenvalue 
problem cannot be recovered by use of this correspondence. Specifically, the limit ya ^ changes the picture 
dramatically in the light mode sector: a pair of the heavy modes from the previous case becomes light and along with 
the already existing pair of light modes, forms a fourfold degeneracy at zero order. On the other hand, the remaining 
pairs of two heavy modes do not undergo a modification and the eigenvectors and eigenvalues stay the same as before, 
with the exception of Smi — > —5m\. 

To summarize, the eigenfrequencies of the pairs of two light modes are 
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and for the pairs of Higgses and longitudinal vectors, we have, respectively, 



= 9l(F,^+Gl) + {3liFoF,+GoG,) 
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where, as before, we identified the longitudinal vector by comparing the frequencies to the mass of A'^^"^ and Af"' 
(ID1[1 at the given order in e. In the above, for later convenience, we defined the dimensionless quantity, 

( :■; rtmS — dm': 1 ( -I- ( tS ' ^ 

AM EE 



[■i 5ml- 5ml) {F^ + Glf 



(GoFq - FoG^f + F^Gl (t^q ~ Eg^ 
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In the degenerate mass limit 6mi — 0, Fq — F, Gq — G, Y,q = E, Eq = E, the eigenfrequencies of the system reduce to 
two copies of the ones of the coupled system in the four field toy model of We will see that this correspondence 
goes even further. 
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The eigenvectors of the hght modes read (see Appendix [E]) 



I I (aj wi + bj W2) cos aj + W2 sin aj 

y/i+a'^+b^. \ (^g^. _|_ Ij . sin aj — W2 cos aj 



U = A, B) 



(D35) 



{aj wi + 6j W2) sin Uj + W2 cos Uj 
{ttj wi + bj W2) cos aj + W2 sin aj 



where the coefficients are 
aA 



bA = bi 



2(f -S) (n-s;, 
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with definitions, 
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The four dimensional vectors wi and W2 in (ID35P are normahzed and orthogonal eigenvectors of matrix A defined in 
(jD2ip . These can be written as combinations of (|D27p through 
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As for the heavy modes, the eigenvectors Vc\i, Vc\2, Vd,i and Vd.2 in (|D26p are also valid for this system. 

The calculation of the matrix T is then straightforward, although the expressions are far from simple. Formally, its 
structure is as follows 
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with a total of ten independent components. We find that, as in the previous system, it is not possible to remove the 
non-adiabatic mixing by rotating the eigenvectors of degenerate states. In general, from the suppression arguments 
made in Section [Till we expect production due to mixings between the light modes and the Higgses. As every physical 
mode has an identical copy, one way to simplify the system is by decoupling the copies into two independent systems. 

We first remove the mixings between the copies of the heavy modes, which are due to the components Fsg and Fyg. 
These vanish by choosing the rotation parameters 



(D40) 



with constant aco and auo- We now simplify the mixings between different modes. For any two sets of distinct 
eigenvalues, there corresponds four components of F matrix. Our goal is to keep only two of these, and remove the 
other two. For the mixing between modes (1, 2) and (5, 6), we choose to solve F15 = which is an algebraic equation, 
with the solution 



tan(a^ — ac) 
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IK 
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Vl + AM2 + 1 
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In a similar way, we remove a pair of components corresponding to the mixings between the modes (3,4) and (5,6), 
by requiring T^q = 0, solved by 

tan(aB ~ ac) ^ \ 9^ 1 ( ^M.^ ) ■ (D42) 




Now, we have fixed all four of the rotation parameters up to integration constants aco and aoo- The logical progression 
of this procedure is to set = and = 0, i.e. to remove the mixing between the two copies of the light modes. 
This, in principle should give two algebraic equations for aco and aoo- However, these equations require time 
dependent solutions, inconsistent with (|D40p . Although we had previously argued that and will not result 
in particle production due to adiabatic behavior of these specific mixings, they will still be responsible for converting 
produced quanta into their twin copy. So our goal of decoupling the two copies cannot be realized in this setting. 
However, relaxing the different mass requirement simplifies the problem considerably. From here on, we will assume 
AM <^ 1 and use the above choices for the rotation parameters and expand the remaining six T matrix components 
in series in AM. At order 0(AM") and 0{e), the remaining components of F are 



Fi2 = 0{AM) , Fi4 = ©(AM) , Fi6 = ©(AM) , F34 = ©(AM) , 
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In this approximation, the modes (1, 3, 5) and (2, 4, 6) decouple from each other and form two copies of exactly the same 
system, which resembles to the coupled system of four field toy model with C/(l) symmetry in By comparing the 
eigenfrequencies, the F matrix components of the two problems are equivalent by the exchange (Fia, F24) — > F23'' 
and (F35, F46) '"° ° i with — > 3e^/2. ^ It should also be noted that the sub-leading terms in the 

eigenmasses of the heavy modes do not match exactly with the ones in the toy model, but being corrections to heavy 
mode mass, they will have a suppressed effect on the final result. However, as mentioned above, in the limit of 
degeneracy, the two problems give exactly the same equations at the given order of e expansion. 

APPENDIX E: DOUBLY DEGENERATE SYSTEM AND UNIQUENESS OF OCCUPATION NUMBERS 

Here, we outline the diagonalization of the matrix form encountered in the last two subsections of Appendix |D] 
The eigenproblem for the 2 N x 2 N matrix of form 

nl + KlK,= (^'^^ ^) , (El) 

where A and B are N x N real matrices with — A and — —B, is equivalent to the eigenproblem of N x N 
Hermitian matrix, 

{A-tB)-iVR + iVi) = X{VR + tVi) , (E2) 

where lambda is an eigenvalue, Vn and Vi are the real and imaginary parts of the corresponding eigenvector. The 
above equation can be written in the form. 



The overall signs of F matrix in llSlI can be recovered by doing a rotation on (1, 2) by an angle (l — sgn[Ff^/Fo — Gq/Gq]) 7r/2 and on 
(3,4) by angle (l - sgn[E^, - %]) 7v/2. 
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The phase freedom in the Hermitian problem, translates to the real case as an SO{2) symmetry on the eigenvectors. 
As there is no way to distinguish between the pairs of the eigenvalues, the eigenvectors can be rotated into one another. 

The problem in the main text is further simplified in the e expansion. In both cases where we have the form (jEl|l . 
we have B — O(e^), and as a consequence, at zero order in e, both vectors Vr and Vj are also eigenvectors of the 
matrix A with same eigenvalue. However, the matrix A has additional degeneracies at zero order, so the eigenvectors 
have another rotational degree of freedom, which can be removed by solving the second order eigenvalue equations. 
In general, the eigenvectors of + KJK^ can be written as 

y 1 I (a wi + 5w2) cosa + W2 sina | y ^ I ^ (awi + bw2) sina + W2 cosa \ 

\/l + + b'^ \ (a wi + 5w2) sina — W2 cosa / ' ^/l + + b'^ \ (awi + 6W2) cosa + W2 sina / ' 

(E4) 

where wi, W2 arc orthonormal eigenvectors of A corresponding to the same eigenvalue, a is the rotation parameter 
arising from the double degeneracy of the fl^ + Kj^K^ matrix. The coefhcients a, b can be determined uniquely 
through the second order eigenvalue problem. In the case where there is no degeneracy in A at zero order, the above 
eigenvectors can be written in a much simpler form 

/ w cosa \ Tr I sina \ 

Vi^ \ . ] , V2= \ , (E5) 

y w sma / \ w cosa / 

where w is the unique eigenvector of A corresponding to the eigenvalue. 

The last point of this section is to show that the occupation numbers are independent of the choice of rotation 
parameters. If matrix ^ diagonalizes fll + KJK4, with {H.^ + KjKi)^^ — diag(a;^), and is arranged such that 
the (2m — l)th and (2m)th columns correspond to eigenvectors of the mth distinct eigenvalue, then ^ = ^TZ also 
diagonalizes it, provided that 

Tlrj = 5rj , (i, j ^ 2to - 1,2m) , 

7?.2)ri-l,2m-l = T^2m,2m—COs9 , 7?.2m-l,2m = ^ 7?.2)ri,2m-l = sin 6* , (E6) 

which rotates the eigenvectors corresponding to the mth distinct eigenvalue, with m — 1, ..N . There are N such 
independent rotations one can apply, under which the matrix F (|A14[) becomes 

f = (^7^ + 7^') . (e7) 

Here and below, an overtilde denotes the quantity calculated using the rotated eigenvectors. Similarly, the transfor- 
mation laws for the matrices / and J pT|) are 

/ = 7^^(/7^ + 7^') , J = n^jn. (es) 

Finally, the Bogolyubov equations will have the same form as in ITU]), 



a , (E9) 



where we defined a = TZ'^a and /3 = TZ^ (]. Of course, as there is arbitrariness in the definition of the eigenstates, the 
occupation numbers of a given rotated state will vary with the rotation angle. However, the physical quantity, which 
is the total occupation number for the mth distinct eigenvalue is invariant of the choice of the rotation parameters 

nm = (r/3^)2m-1.2m-l + (/3'/3^)2rn.2m = (7^^/3*/3^7^)2™-l,2m-l + {Tl^ f3'^Tlhm,2m = n„, . (ElO) 
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